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The routing capability is a requisite in quantum network. Although the quantum routing of signals has been 
investigated in various systems both in theory and experiment, the general form of quantum routing with 
many output terminals still needs to be explored. Here we propose a scheme to achieve the multi-channel 
quantum routing of the single photons in a waveguide-emitter system. The channels are composed by the 
waveguides and are connected by intermediate two-level emitters. By adjusting the intermediate emitters, 
the output channels of the input single photons can be controlled. This is demonstrated in the cases of one 
output channel, two output channels and the generic N output channels. The results show that the 
multi-channel quantum routing of single photons can be well achieved in the proposed system. This offers a 
scheme for the experimental realization of general quantum routing of single photons. 

Quantum routing^ of information from one sender to many receivers is an essential function for the 
quantum network^^. Single photons are suitable candidates as the carrier of quantum information due to 
the fact that they propagate fast and interact rarely with the environment. The field-matter interaction 
can be used to manipulate the photons. By the way, the light-light interaction by the platform of linear optics^'^'^^'^^ 
can also be used to manipulate the photons. The photon transport in a one dimensional (ID) waveguide has been 
studied extensively both in theory^^"^^ and in experiments^^"^^ because the strong coupling of the waveguide- 
emitter can be achieved. In the waveguide- emitter system, the waveguides act as the channels and the emitters as 
the nods of the quantum network. Based on these advantages, the quantum routing of photons in the waveguide- 
emitter system is promising. Recently, two output channel quantum routing of single photons is studied the- 
oretically in a waveguide- emitter system^'^. The input single photons can be redirected into either of the two 
output channels with a maximal probability of unity and no more than 1/2, respectively. It is interesting if the 
photon can be redirected into either of the output channels with an extremely high probability. Moreover, a 
general form of many output channel quantum routing of the single photons in the waveguide QED system will be 
of considerable interest. 

For these purposes, we propose a novel scheme for the quantum routing of single photons from one input 
waveguide into N output waveguides. In our scheme, the ith output waveguide is connected with the input 
waveguide by an intermediate two-level system (TLS). After scattering, the single photon injected into the input 
channel is redirected into other channels. For the single output channel quantum routing, the quantum inter- 
ferences redirect the input photon into the output channel completely when the intermediate TLS resonantly 
interacts with the input and output channels with the same strength. The single output channel routing properties 
can be modified when an additional TLS is coupled to the input channel. In the two output channel case, the 
photon can be redirected into either of the output channels with an approximate unity probability. The photon 
can also be redirected completely into both the output channels with various probabilities. In the generic N output 
channel case, the quantum interferences prevent the photon being redirected into the other channels except the 
input channel when all the TLSs resonantly interact with the channels with equal strengths for a large value of N. 
They also can completely prevent the photon being directed back into the input channel for suitable parameters 
for any value of N. The photon can be redirected into a desired channel with an approximate unity probability. 

The schematic diagram of the system under consideration is shown in Fig. 1. The input channel is a Sagnac 
interferometer^^"^^ composed by a semi-infinite waveguide, a 50 : 50 beam splitter and a waveguide loop. The 
input channel connects the ith infinite waveguide by an intermediate TLS at the middle point of the waveguide 
loop. For simplicity, we label the interaction position x = O.lt is convenient to bring in the even and odd operators 

as Uek = {cLrk + ciik) and Uok = {cirk — ciik)y with (ciik) being the annihilation operator for the right (left)- 
V 2 V 2 
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Figure 1 | Schematic diagram of the multi-channel quantum routing of 
the single photons. AT ID waveguides are connected with the input channel 
by N intermediate TLSs. 

moving mode with the frequency Vg\k\ in the waveguide. Hereafter 
we will take the photon group velocity Vg = 1. The clockwise- and 
counterclockwise-moving photons in the waveguide loop can be 
treated as the left- and right- moving photons. If we inject a photon 
into the input waveguide, an even mode quasi particle can be pro- 
duced at the middle point of the waveguide loop when the phases of 
the clockwise- and counter clockwise-moving photons are equal. In 
the even and odd picture, the odd mode only contributes to the free 
energy part in the Hamiltonian. Therefore, it is enough to study the 
dynamics of the even mode. For simplicity, we omit the superscript e 
of the even operator a^k- The even part of the Hamiltonian can be 
written as {h = l) 

pco r pco 

H= dkka\ak+ | | dkka\^ai^k + 0JiG^ a^' 

J — GO y L J — GO 

(1) 

/•GO 

J — 00 

with a\ being the even mode creation operator in the input (/th 

output) waveguide, and oji being the /th TLS transition frequency. 
We have taken the energies of the TLS ground states zero. The terms 
of the second line in Hamiltonian (1) represent the interaction of the 

TLSs with the waveguides, -yz and are coupling strengths of the 
v 2 v 2 

/th TLS to the input waveguide and the /th output waveguide, 
respectively. The coupling strengths can be written as 

-(C = ± ), with y] being the decay rate from the /th TLS to 

the waveguides. The TLS can be a real two-level atom or an artificial 
two-level system. In the single-photon case, the TLS can also be a 
cavity or a cavity- atom dressed system. The transition frequency of 
the TLS can be tuned by Stark shift and the coupling strength 
between the waveguide and the TLS can be varied by changing the 
distance between them^°'^\ Especially, for a waveguide -atom dressed 
system, the effective coupling strength of the dressed TLS to the 
waveguide can be adjusted by the extra laser^^. We have made two 
approximations: one is to extend the frequency integration to ±0°, 
the other is to assume that the coupling strengths are independent of 



the frequencies, which is equivalent to the Markovian approxi- 
mation. These approximations are valid since we will focus on the 
pulse with a narrow frequency width around the carrier frequency. 
The arbitrary state of the system in the single -excitation subspace 

hastheformof |^^A;afV)4+ \^^'^tk^M,k + 

{iW ] l<^>' with ^ (0, (0 and ^^^^ (t) being the probabil- 
ity amplitudes. The state |(/>) represents that all the waveguides and 
TLSs contain no excitation. We assume that, initially, a photon pre- 
pared in a wave packet with a Lorenzian spectrum is injected into the 
input waveguide, while the TLSs and the output waveguides contain 

/ejn 



no excitation, i.e. oc?^^ (0) = 



e and are the spectral width and the center frequency of the input 
wave packet, respectively. When the width e^O, the frequency of the 
carrier is equal to the central frequency. This corresponds to the 
monochromatic limit. 

Results 

Single output channel. When the input channel is connected with 
only one output channel by a TLS, the probability amplitudes are 



1^ I ^ 

2^2 



and a, 



Obviously, only when = 1 and = 0, the input photon is 

redirected into the output channel completely. This is due to the 
quantum interferences. When the detuning is large enough, the 
photon will be back into the input channel with an approximate unity 

maxjy^ ,y r | 

probability. Here we bring in a parameter = ^— — ^ to 

min{y+,y,- } 

measure the difference between the coupling strengths of the /th 
TLS to the input and the /th waveguides. When ^i^l, we can get 

^ ~1, which is similar to the large detuning case. Although the 

decay rate of the TLS to the input channel is much smaller than the 
decay rate to the output channel, the photon is prevented being 
redirected into the output channel. The condition f/i ^ oo 
corresponds to the limit that one of the two waveguides is 
decoupled to the TLS. In this case, the input and output channels 
are not connected, and the photon will be back into the input channel 
completely. To show the details of the quantum routing of the single 

phonon in theN = 1 case, we plot the probability (x[ 



2^2 



, with di 



Wi 



(1) 



against the 
^ IS not 

I 1 2 II 

plotted here because a[^M = 1 — k^/^ . The probabilities of the 



detuning and decay rates in Fig. 2. The probability 

i2 

I --^~ru 




Figure 2 | Probability of the single photon in the 1st output channel in the 

I (1)1^ 

long-time limit ^ against the detuning and decay rate when N = 1. 
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photon in the input and output channels can be controlled by 
adjusting the detuning and coupling strengths. The single photon 
can be directed into either of the two channels completely or directed 
into the two channels with a desired probability. The N = 1 case is 
analogous to a waveguide coupled to a A-type three-level quantum 
emitter^^'^^''*^ 

Two output channels. When the input channel is connected with 
two output channels by two TLSs, the probability amplitudes are 



obtained as a 



iS2 + ■ 



( 



)■ 



4 



Because 



and 



the 



(2) (2) 

expression of a^^ is similar to the study of the properties 

I (2) P I (2) P 

^ and is enough. The large detuning between the 1st 

TLS and the input photon prevents the photon being redirected 
into the 1st channel. Especially, when ^2^0 and y+ =y~, it also 
prevents the photon being directed into the input channel. As a 
result, for an extremely large value of ^i, the photon will be 
redirected into the 2nd channel with an approximate unity 
probability when ^2 — 0 and y2 • And when both the 
detunings and 62 are large enough, the photon will be directed 
into the input channel with an approximate unity probability. 
Besides, the photon distribution can be influenced significantly by 
the coupling strengths. For example, when ^2 — 0 and y2 ^ the 
large value of y^ prevents the photon being directed into both the 
input and 1st channels. When r]i = f]2 = ^ and ^1 = ^2 ^ 0^ the input 
photon will be directed into the input channel with a small 
probability and be redirected into the other two channels averagely 
with a large probability. To see the details, we plot the probabilities 

I (2) P I (2) P 

^ and against the detuning and coupling strength of the 

2nd TLS when ^1 = 0 and y^ = y^ = y^ in Fig. 3(a) and 3(b). It is 
interesting that when ^2 and y^ is small enough, the photon will be 
back into the input channel with an extremely high probability. This 
will be studied in detail below. 



If we do not wish the input photon to be back into the input 
channel, it is easy to choose the appropriate parameters which satisfy 

(2) 

— 0. Here we take two simple cases to investigate it. One case is 
when y^ = y^ and y2 =72 ^ the condition 4^i ^2 + Yi 72 =^ ^an be 
satisfied for various values of the detunings. The other case is when ^1 
= 62 = Oy the relation {y^ ~yr) (y^ "^2") ^^1^2" easily 
satisfied. Especially, in the latter case, when yi = 2y^ and y2 = 2y2 , 
the input single photon is redirected into the two output channels 

with equal probability - . We plot the probabilities of the photon 

being redirected into each of the two output channels in the two 
cases in Fig. 3(c) and 3(d). It shows that the input photon can be 
completely redirected into the output channels with various 
probabilities. 

It is necessary to study a special case of AT = 2, that is, the 2nd TLS 
is decoupled to the 2nd output channel, i.e. y2 = 0. This can be 
understood that an additional TLS is coupled to the input channel 
in the AT = 1 case. The additional TLS will modify the system beha- 
vior. For example, when r]i = 1 and ^1 = 0, we find 



(2) 



-/4^2yi 



and a[^^ = - 



Hence, when 



4/^2?! +yi 72 " ^ 4/^271 +yi 72 ' 
the additional TLS resonantly interacts with the input photon, i.e. 
^2 = 0, the photon is directed into the input channel compared with 
the AT = 1 case, in which the photon is redirected into the 1st output 
channel. This can be seen in Fig. 3(a) and 3(b). When 62 is large 
enough, the input photon will be almost completely redirected, 
mapped to the AT = 1 case. By adjusting the additional TLS, the 
photon distribution in the input and output channels can be con- 
trolled. This provides more control to the single output channel case. 

N output channels. We proceed to study the general case that the 
input channel connects with AT output channels by AT TLSs. Let's first 
consider the simplest case that all of the N output channels are 
identical. That is to say, all the TLS are identical, all the decay rates 
to the input channel are identical, and all the decay rates to the output 
channels are identical. We label 3i = S,y^ =y^ , and y^ = y~ . The 



probability amplitudes are obtained as a[^^ = 



, and 



2/^ + y+ -Ny~ 
2id + y+ +Ny- 



^V7~7^ 
2i3 + y+ +Ny- 

It is noted that when all the TLSs interact 





Figure 3 | Probabilities | afl | , | afl | and | a[^^ | against the detunings and coupling strengthes when N = 2. (a) and (b) are the probabilities | afl | and 

I f2) I ^ I f2) I ^ I f2) I ^ 

a[ ^ when ^1 = 0 and = y^ = y^ , respectively, (c) and (d) are the probabilities W^A WiA '^'^ the probabilities against the detuning ^1 when 
y^ = yj~ , y2^ = y~ , and 43 162 + y^ y2 =0.(d) is the probabilities against the coupling constant y^ when di = 62 = 0, and (y^ — y^ ) (y2 ~ ) = Yi" 72 - 
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resonantly with the input photon and the coupUng strengths satisfy 
y+ = Ny~y the interferences prevent the input photon being directed 
into the input channel and redirect the input photon into each of the 

N output channels with equal probability — . When all the coupling 

strengths are equal and all the detunings are zero, the probabilities 



are obtained as 



and 



As 



1-N\ 

i+ivy 

redirect the input 



discussed above, the quantum interferences 
photon completely from the input channel into the other channel 
when N = 1. As the number of the output channels increases, the 
probability of the photon back into the input channel increases. 
When N = 3, the input photon is distributed in each of the four 
channels, including the input and output channels, with equal 

probability -. When the number of the output channels is large 

enough, the quantum interferences direct the input photon back 
into the input channel almost completely. 

We now study a simple case to illustrate how to redirect the input 
photon into the desired channel. We assume that AT — 1 of the N 
channels are identical except the mth channel, we label 3i^^ = d\ 
and yf-^^=y~'. The probability amplitudes are obtained 



(N-l)y- 



and 



(N-l)y- 



When 



{iS„+f){i6'+' 
the decay rates y ~ = are much larger than other parameters. 

The photon will be redirected into 



2 ' 4 

the mth channel almost completely when 2y^'^{N—1)y~'. The 
limit of this condition is that the iV — 1 TLSs are decoupled to the 
input channel. Besides, when all the detunings are large enough, the 
photon will be directed into the input channel with an approximate 
unity probability. When Y~ =Y^,Sm = 0, and 3' is large enough, the 
photon will be redirected into the mth channel with an approximate 
unity probability. Similar to the N = 2 case, the photon can be 
redirected completely into the output channels where the photon 
probabilities can be controlled. As a consequence, the N output 
channel quantum routing of the single photons can be achieved. 

The single photons are considered as the ideal carrier of quantum 
information. The quantum information can be encoded by the rou- 
ted photons. We label two degenerate photon modes as ai and ^2. For 
example, the two degenerate modes can be horizontally and vertically 
polarized modes or can be right and left circularly polarized modes. 
To achieve the quantum routing of the quantum information, we 
bring in N extra intermediate TLSs to make sure that each output 
port is connected with the input port by two intermediate TLSs. The 
distance between the two intermediate TLSs is much smaller than the 
wavelength of the photon. One of the two intermediate TLSs only 
interacts with ai-mode photon and the other only interacts with 
mode photon. In this case, the Hamiltonian can be written as 



H 



dkka 



t 



■ r 



^mi+ J-oo ^K^mi^rni,k+g~iamk)o^i + h.c]]. To study the quantum 
routing of the quantum information, we take the input state as 

\^>)=[A^f\Q)al + Baf\Q)a\^\i>). with |A|^ + |B|^ = 1. 

Obviously, when^^l" ^^it' ^Sii ^iv the input state 

1^) can be routed to the desired channel with desired probability. In 
this case, the information remains unchanged after being routed. In 
our scheme, the control information is stored in the coupling 
strength and the TLSs. It will be interesting to use quantum 



information to control the routing. This might be achieved by pre- 
paring connections between the control quantum states and the 
target quantum states. 

The authors in Ref 9 have raised a set of five requirements 
imposed on a full quantum router. It will be interesting when a router 
meets all the requirements. As discussed above. Our scheme will 
meet all those requirements if the control information can also be 
encoded into quantum states. This full quantum routing scheme in 
various systems is an interesting question which need be further 
studied. 

Discussion 

We have investigated the quantum routing of the single photons 
from one channel to multiple output channels in the proposed sys- 
tem. We have shown that the input photon can be redirected into any 
of the channels with a high efficiency in our scheme. Alternatively, by 
cascaded combining many two output channel routers demonstrated 
in the presented literatures, the multiple output channel quantum 
routing can be achieved. We note that our general N output channel 
quantum routing shows essential advantages compared with the 
cascaded combined one. In the cascaded combined scheme, the effi- 
ciency of the / output channel quantum routing is limited by the 
efficiency of the / — 1 output channel quantum routing. However, 
our scheme is immune to these drawbacks. It is necessary to point out 
that our scheme is equivalent to that the input waveguide loop is 
coupled to N output waveguide by an intermediate N + 1 -level 
system at the middle point of the waveguide loop. Our scheme sheds 
light on the experimental realization of quantum routing. 

Methods 

The continuum form of Hamiltonian (1) can be derived from its discrete form by 
following the technique in Ref. 45, 46. The free Hamiltonian of the photon in the ID 
cavity of length L has the form of Hq = <jJka\ak, and the interaction Hamiltonian 
of the photon with the TLS has the form of Hi = g^l^~ + h.c. The sum can be 
converted into an integral by — dk. As done in Ref. 45, we can bring in the 



operator a 



-at and obtain the Hamiltonian (1) with continuum form. 



The Shrodinger equation governed by Hamiltonian (1) gives a set of differential 
equations of probability amplitudes. By performing the Laplace transformation, the 
set of differential equations under the initial condition are transformed as 



/n/e 



{s-\-i(i)k){(i)k — zu-\-ie) 



-JSi_dN) 



s + icoi^k 



(2) 



j=l..N 



s+e+iw 



We proceed to obtain the expression of (s). For simplicity, we label 



a,- = s + ioji H — — , hi - 
2 

can be written as 



y- ^Tip"- v/Ti/e 

I and Ci = — ■ — • "^^^ ^^^^ equation in equations (2) 



(3) 



(N) 



j=l..N 

From equation (3), we can find ^ bjl^j^\^) = ~ n ~ 

J2j=i..N h^f^ (^) - Ef = i..]v to obtain the expression of ^.^^ ^ hj^] 
Consequently, the expression of (s) can be obtained as 



5+e+ iw 



(4) 



From the equation (4) and equations (2), the expressions of a[^^ (s) and (if^ (s) can be 
easily found. 
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In the long-time limit, the amplitudes can be found by performing the inverse 
Laplace transformation of a[^^ (s), a|^^ (s) and pf^'^ (5). Then we prove that the 

amplitude (t) ~0 when t^^.lt means that the TLSs are in their ground states in 
the long-time limit due to the interaction with the continuum of modes. Obviously, 

n^.i..„(5+,-co^+y+) + 5^,^^^ [r-Uf^j^s+mf+Y;)] = ^J^Js-a,), 

with ai being the root of the equation 

llj=i..N (5 + icnj + y/ ) + Y.j=i..N W ^7 + ^'w/ +?/)]= 0. We lable a, = x, 
+ iyi, with Xi and)/, being real numbers. We will show that the value of Xi can only be 
smaller than zero. 
If Xj > 0, we can get 

^j=i..N (at + iojj + y + yY.j=i..N W (^'- + + ) ] 

= 11^= {irrij + rij) + J2j=Q..N [yf {^^f + )] = 0' ^ith m^- and being 

real numbers and rij > 0. We label imj + rij = Rje'^^ and H^- = i.N(/mj + = Re''^, and 

get 



;=1..N 



j=l..N 



3n 



Obviously, it can not be zero because -> 9j>0 or 2n>9j> and 
|y^~ , R, Rj I > 0. Therefore, the value of Xj must be smaller than 0. As a consequence, 
wecanfindjSS^^(0 =0. 

lim 00 

Similarly, the probability amplitudes of the photon in each channel in the long- 
time limit are found as (0 = a[^^af ^ (O)^-'''^^ and af ^ (t) = ^af ^ {0)e-^'\ 

limt^cx) lim 00 

with 



}=i..n\ ^ 2y .^/-.^ 



^ n ( iSf + ^ 



n 

;=1..N 




- E 

j=l..N 


'in 

. f*i 






n 

j=l..N 




+ E 

; = 1..N 


in 







(5) 



We have taken the conversation of energy condition cOi^k = (^h After scattering, the 
probabilities of the single photon in each channel are a|^^ and a[^^ . 
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